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ABSTRACT 


The  propagation  of  infinite  trains  of  symmetric  harmonic  waves 
traveling  in  infinitely  long,  right  circular  cylindrical  shells  is  investigated 
on  the  basis  of  the  three-dimensional  theory  of  elasticity.  The  shells  are 
assumed  made  of  three  concentric,  transversely  isotropic  cylinders,  each 
of  different  materials,  bonded  perfectly  at  their  interfaces.  The  frequency 
equation  is  established  by  representing  the  displacement  field  in  each 
cylinder  in  terms  of  potential  functions  and  satisfying  the  Navier  equations 
of  motion  and  the  boundary  and  interface  conditions  of  the  cylinder. 

The  frequency  equation  has  been  programmed  for  numerical  evaluation 
on  an  IBM  7044/7094  DCS  computer,  and  the  influence  of  the  mechanical 
properties  of  the  layers  on  the  frequencies  of  the  first  few  modes  is  investigated. 
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Cylindrical  coordinate 
Non-dimensionaly zed  radia1  coordinate 
Radii  of  the  shell  layers  (see  fig.  1) 
Nondimentionalyzed  Radii  of  the  shell  layers 
Thickness  of  i**1  layer 

displacement  of  components  in  the  r  and  z 
respectively 

Component  of  strain,  defined  in  eq.  f  l]  . 
Component  of  stree  defined  in  Eq.  f  l]  . 
density 

Elastic  constants 
Nondimensional  elastic  constants 
Wavenumber  in  the  axial  direction 
non-dimensionalized  frequency 
circular  frequency 
non-dimensionalized  wavenumber 

Time 

"veloc  ities” 

nondimensional  velocity  ratios 
radial  wave  numbers  (see  equation  [8]) 
see  equations  [19]  and  [20] 


zero  order  Bessel  function  of  the  first 
or  second  kind  respectively 
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INTRODUCTION 


The  frequency  equation  for  harmonic  waves  traveling  in  traction-free 
infinitely  long,  isotropic  circular  cylindrical  shells  has  been  established 
on  the  basis  of  the  three-dimensional  theory  of  elasticity  and  has  been 
evaluated  numerically  by  Gazis  (4)*and  Greenspon  (5).  More  recently, 

Mirsky  investigated  the  propagation  of  harmonic  waves  in  circular  cylindrical 
shells  made  of  transversely  isotropic  and  of  orthotropic  matericals  (9)  (10). 

The  increasing  demand  for  structural  components  of  aerospace  vehicle's 
having  a  high  strength  to  weight  ratio  and  being  capable  of  withstanding 
high  temperatures,  has  resulted  in  extensive  use  of  multi-layered  shells 
and  in  considerable  interest  in  the  propagation  of  harmonic  waves  in  such 
shells.  ArmcnSkas  (1)  (l)  presented  a  unified  treatment,  on  the  basis  of 
the  theory  of  elasticity,  for  harmonic  waves  of  an  arbitrary  number  of 
cicumferential  nodes  traveling  in  two  layered  isotropic  shells.  Keck 
and  Armenikas  (7)  investigated  the  propagation  of  harmonic  axisymmetric 
waves  in  sandwich  isotropic  shells.  Moreover,  a  number  of  approximate 
theories  for  two  and  three  layered  shells  were  established  (3)  (6)  (11). 

In  this  investigation,  the  frequency  equation  for  propagation  of 
trains  of  axisymmetric  nontorsional  harmonic  waves  in  infinitely  long 
shells,  made  of  three  concentric  cylinders  of  different  transversely 


2 


isotropic  materials,  is  derived  on  the  basis  of  the  linear  theory  of 
elasticity.  It  is  shown  that  as  in  the  case  of  isotropic  shells,  for 
waves  having  infinite  axial  wave  length,  the  frequency  equation 
degenerates  into  two  independent  equations  for  uncoupled  longitudinal 
shear  and  uncoupled  radial  motion. 

The  frequency  equation  has  been  programmed  for  numerical 
evaluation  on  an  IBM  7044/7094  DCS  computer,  and  the  influence  of  the 
mechanical  properties  on  the  frequencies  of  the  first  few  modes 
is  investigated. 
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SOLUTIONS  OF  THE  EQUATIONS  OF  MOTION 


In  the  ensuing  derivations,  Che  following  notation  for  the  components 
of  stress  and  stain  will  be  used 

Tll“Tl’  T22“V  T33_T3*  T23“V  X3l“V  T12“V  ^ 

en“er  e22“e2*  e33“e3*  2e23“e4*  2e3l“e5*  2e12"e6  * 

In  terms  of  this  notation,  Hooke's  law  for  a  general  anistropic  body 


may  be  written  as 


'1  ■  ‘ij'j  <*•  i  ■  !•  2 . «• 


For  a  transversely  isotropic  body,  in  particular,  these  relations 


reduce  to 


c66e6  ’ 


where 


Cliei  +  C12e2  +  C13e3* 
c12ei  +  clle2  +  C13e3’ 
C13C1  +  C13e2  +  C33e3’ 
W  T5  ‘  c44e5’  T6 


C66  "  !j(cll  -  C12)- 


The  assumption  that  the  strain  energy  density  is  a  positive  definite 

quadratic  function  of  the  components  of  strain  imposes  the  following 

restrictions  on  the  elastic  constants 

c  >  0,  c  >  0,  c  >0, 

11  33  44 

C11  ‘  c12  *  °*  C11  C33  "  C13  ”  °*  ^ 


C11  C33  +  C12  C33  ”  2C13  ”  °* 
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Equations  [3]  reduce  to  those  for  an  isotropic  body  by  empl tying  the 
following  relations  between  the  elastic  constants 

C33  "  cll*  c12  "  C13'  c44  “  c66* 

The  components  of  stress  may  be  obtained  in  terms  of  the  components 
of  displacement  by  substituting  the  strain-di  -placement  relations  into  the 
const i tuat ive  equations  [3].  These,  in  turn,  may  be  substituted  into  the 
stress  equations  of  motion  to  obtain  the  following  displacement  equations 
of  motion 

cU(ur,rr  +  ur.r/r  ’  “r/r2)  +  c66  ur,ee/r2  +  c44ur,«  + 

2 

(c.  +  c,„)u  / r  -  (c,,  +  c,,)u  /r  +  (c  +■  c  )u  ■  pu  , 

6b  12'  e,re  66  ll'  6,6  13  44  z.rz  w  r 


(c!2  +  C66>Ur.re/r  +  <C11  +  C66)ur.8/r2  +  c66<Vr  +  “e.r/r  * 


+  C,,  Un  aJT  +  C//UQ  +  <C1  1  +  CA/)U,  *JX  “  PUft» 

11  6,66  44  0,zz  13  44  z,6z  o 


(c._  +  c..)  (u  +  u  /r  +  u  /r)  +  c  u  + 
13  44  r,rz  r,z  6,rz  13  z,zz 


c  (u  +  u  /r  +  u  aJr  )  -  pu  . 
44  z, rr  z, r  z,66  z 


Here  p  is  mass  density  per  unit  volume;  subscripts  preceeded  by  a 
comma  denote  differentiation  with  respect  to  the  space  coordinates.  The 
dot  indicates  differentiation  with  respect  to  time.  It  can  be  shown  (8)  that 
for  axisymmctric  motion,  Eqs.  [5]  are  satisfied  by  a  displacement  field 
of  the  following  form 


uf  -  (c13  +  c44)*(r)  r  cos (kz  -  *t), 


6  J 
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i  2  ,2 

u  ■  —  [  e , ,  V  $(r)  +  -  c  )  $(r)]  sin  (ka  -  at),  [6  cont] 

z  k  ii  1  44 


where  the  potential  function  $(r)  must  satisfy  the  equation 


+  p2)  (V*  +  q2)  4>(r>  -  0. 


Here 


;2  -  d  +id 

l 


and  u  and  k  are  the  frequency  and  the  axial  wave  number  of  tne  wave, 
respectively.  The  radial  wave  numbers  p  and  q  are  given  by 


2C11C44 


[A  i  /A  -  B  ]  , 


wnere 


A  -  cn(.  -  c„)  +  C//  (m  -  cLt)  +  (c^  +  c^)  , 


"  "11'“  "33'  T  44'  "447 

B  "  4C11  C44(n  -  C33>  (m  “  c44}  • 


2..  2 

a  ■  pu  /k  . 

2 

The  plus  sign  in  equation  [7]  refers  to  p  .  Notice,  when  the  radicand 

2 

A  -  B  is  negative,  that  the  radial  wave  numbers  p  and  q  become  complex.  Thus 
for  a  certain  range  of  the  values  of  the  elastic  constants  there  is  a  range 
of  values  of  real  m  for  which  the  radial  wave  numbers  p  and  q  become  complex. 
In  the  case  of  isotropic  elastic  shells,  the  radial  wave  numbers  do  not 
assume  complex  values  for  any  real  value  of  m.  Thus,  it  seems  appropriate 
to  classify  transversely  isotropic  materials  as  (a)  less  anisotropic  if 
their  mechanical  properties  are  such  that  p  and  q  do  not  become  complex  for 
any  value  of  m,  (b)  more  anisotropic  if  their  mechanical  properties  are  such 
that  p  and  q  may  become  complex  for  a  certain  range  of  m. 
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The  operators  in  equation  [7]  are  Bessel  operators  and,  conse¬ 
quently,  their  solution  is  given  in  terms  of  the  rero  order  Bessel  functions 
of  the  first  and  second  kind  with  real,  complex  or  imaginary  argument3, 
depending  on  whether  p  and  q  are  real,  complex  or  imaginary.  Therefore 
in  order  to  specify  the  real  solutions  for  $,  and  the  displacement  and  stress 
fields,  it  is  necessary  to  establish  the  range  of  the  material  properties 
and  the  range  of  values  of  the  wave  parameters  (positive  values  of  u  and  £) 


for  which  p,  q  are  real,  imaginary  or  complex. 

2 

The  radicand  A  -  B  of  equation  [8]  vanishes  if  m  satisfies  the  equation 


nJ 


“  %4  )(C11  "  ‘W  "  (C11  +  C44)(C13  +  C4h)2 


2(C13  +  C44) 


(cll  "  c44} 


(cll  _ 


[10] 


/c,,  c 


11  44 


t<C,,  +  C44>  -  <cl,  *  c,./><c- 


-  c, 


13  *44'  ''11  '44' ''33  '44 


)] 


For  ordinary  engineering  materials,  it  may  be  assumed  that  c  >  c  .  On 

33  44 

2 

this  basis  it  can  be  shown  that  the  radicand  A  -  B  cannot  vanish  for  values 

2 

of  m  greater  than  c^.  Thus,  for  a  given  material,  the  radicand  A  -  B 
becomes  negative  for  values  of  m  satisfying 


m  <  m  <  in  (  c 
1  2  44 

0  ^  m  <  m-  <  c 

2  ‘  44 


if  m^  >  0, 


if  m^  <  0. 


If  the  inequality,  A  -  B  <  0,  is  solved  for  c^,  the  following  inequality 


results 


lc13  +  c44>  '  /c44<c44  '  m)r 
<C33  -  m> 


<  c  < 
'  11‘ 


[(c13  +  c44>  +  /c44<c44  -  m,’)' 
<c33  ‘  ") 


.  [12] 
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The  cross-hatched  region  in  Fig.  1  represents  the  locus  of  the  values 
of  m  and  Cjj  satisfying  this  inequality.  Note,  that  when  m  =  0,  inequality 
[12]  reduces  to 


2 

13 

<  r 

,  ,c13  +  2c44»2 

C11 

33 

c33 

i  m 

•*r 

•>r 

u 

II 

inequality  [  1 2  ] 

(c13 

+  C44>2 

[13] 


[  14] 


C33“C44 


The  value  of  m  corresponding  to  the  maximum  value  of  Cjj,  for 
2 

which  the  radicant  A  -  B  vanishes,  can  be  obtained  by  setting  the  derivative 
of  rn  with  respect  to  m  (see  fig.  1)  equal  to  zero. 


d(C-n)  d 

dm  dm 


Ac13tc44ttc44(c44-m|P]  ]  ,  „ 
\  - m  / 


[15] 


This  results  in 


=  c44  1 


I  for  Tj  j  max 


(c33  ’  C44) 


(c13  +  c44) 


(T, , )  =  c  . .  + 

Umax  44 


(c13  +  c44> 
c33  ’  C44 


[16] 


From  equation  [  10  ]  it  can  be  seen  that  the  maximum  value  of  TTj  j 

is  also  obtained  when  m  =  m^  =  m^.  From  the  aforegoing  discussion  we 

may  conclude  that  if  the  elastic  constants  of  a  material  satisfy  the  relation 

(c._+  c..)2 

c,,*c44+  11  44  .  [14] 

c33  ' c44 
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the  radicant  A  -  B  is  positive  for  all  values  of  m  and,  consequently,  the 

radial  wave  numbers  p  and  q  do  not  become  complex  for  any  values  of  m. 

In  this  case,  the  material  has  been  classified  as  less  anisotropic.  If, 

however,  the  elastic  constants  of  a  material  do  not  satisfy  inequality  f  17  ] , 

2 

then  for  a  certain  range  of  values  of  m  the  radicant  A  -  B  will  be  negative 

and  p  and  q  will  be  complex.  In  this  case,  the  material  has  been  classified 

as  more  anisotropic.  Referring  to  equation  [8],  it  can  be  seen  that  p^ 

2 

is  the  conjugate  of  q  .  Consequently,  two  pairs  of  complex  values  of  p 
and  q  are  found.  It  can  be  shown  that  any  two  of  the  four  roots  satisfy 
the  requirements  of  the  equations  of  motion.  The  value  of  q  used  in  this 
analysis  is  the  negative  complex  conjugate  of  p. 

From  Eqs.  f  9],  it  can  be  deduced  that  for  m  >  c^  the  radial  wave 
numbers  p  and  q  are  real,  while  for  c^  >  m  >  c^,  p  is  real  and  q  is 
imaginary.  Finally,  for  c^  >  m,  the  non-complex  values  of  p  and  q  are 
both  real  if  A  >  0  or  both  imaginary  if  A  <  0.  The  parameter  A  vanishes  if 
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the  parameter  A  is  negative  for  all  values  of  m  smaller  than  c  .  For 

44 

these  materials,  for  >  m,  the  non-complex  values  of  p  and  q  are 

imaginary.  For  more  anisotropic  materials  with  elastic  constants 

2 

satisfying  the  relation  c,,  <  (c  +  c,.)V(c._  -  c,,),  the  parameter  A 

11  13  44  33  44 

is  positive  for  c  >  m  >  m^.  Therefore,  in  this  range,  p  and  q  are  both 
real.  On  the  basis  of  the  foregoing  discussion,  the  radial  wave  numbers 
p  and  q  may  be  written  in  the  following  form 

p  -  (e^*5  ko-  q  -  k/3.  [  19 J 

For  less  anisotropic  materials  and  for  more  anisotropic  materials,  if 


q  -  U2r  k/3 . 


m  <  or  m  >  m^,  the  parameters  a  and  6  are  given  by 


1-— i 
4 


A  t  /  A  -  B 


11  C44 


and  alternatively,  for  more  anisotropic  materials  with  m^  <  m  <  m^  (here  m 
replaced  by  0  if  <  0)  a  and  0  are  taken  as 


with 


T**4  (±  cos  0  +  i  sin  0)  , 


?  -  — <»■  -  cJ3)  -i-  <m  -  c  >  , 


0  -  ^  tan"1 


'B  -  A 
A 


For  materials  with  c^  >  (c^  +  c  )  /  (c^  “  »  t{ie  si8n  factors 

(i-1,2)  are  assigned  the  values 


r  1 

if 

m  >  c 

44 

-1 

4 

if 

c  >  m  a  m 

44  2 

1 

if 

in  >  m  >  ra, 

2  1 

.-1 

if 

m^  >.  m  *  0, 
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Consider  a  shell  made  of  three  concentric  cylinders  of  different  materials 
perfectly  bonded  at  the  interfaces.  The  solution  obtained  in  the  previous 
section  may  be  applied  to  each  layer  of  the  shell.  The  material  properties 
of  the  layers  will  be  identified  with  the  superscripts  "(1)",  "(2)",  and 
"(3)"  for  the  inner,  middle  and  outer  layers,  respectively  (see  Fig.  2). 
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For  convenience,  the  following  non-dimensionalized  parameters  are 


introduced. 


r,  a,  b,  c,  d  »  (r,  a,  b,  c,  d)  ~  .  ;  »  kh^/n  , 

h(2) 


n  -  uh(2)  v  P(2)  h  ‘v'c44(2)  .  uj  “  uj  (2)  (J  °  r*  2)* 


,(i)  cn 
1  "  r  (2) 
44 


.(i)  c12 

2  “  7T 7) 

44 


,(D  ci 


|(i).  _A4_ 


4  c  (2) 
44 


d(i)  - 

d5  (27 

C44 


(i  =  1,2,3)  [26 


(i)  _  „  (i) 


d6^  — 


, 


The  ratios  [c  p^]  and  [c/^/  p^]  are  referred 

44  33 


to  as 


"velocities"  and  are  denoted  by  v/l^and  vrf  ,  respectively.  In  non-dimension¬ 
alized  form,  these  velocities  will  be  taken  as 


a<°-  > 

s 


an,  .(i)  _  vd 

and  b  ■  —  (2)  * 


(i  -  1,2,3).  [27] 


The  radial  wave  numbers  may  now  be  re-written  in  terms  of  tiie 
non-dimensionalized  parameters  as 


-(i)  /“ID  (i) 

P  “  /  4 CO 

-(i)  r~\ T)  (i) 

a  -  /  e„  '.B 


where,  for  less  anisotropic  materials  and  for  more  anisotropic 
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materials  if  m  <  m^  or  ra  >  the  parameters  a^and  i^^are  given  by 


.  —  *(i>  .  A«)2  .  ,(!)' 


a)  rr 


with 


- 


Ai)  . 

B  *  4 


2  \  H(0  2  \  H(0  /h(1) 

£  . ,  u  ‘>2 .  .  A  +  it- 

„<*>\*  )  j  dHd<*> 


(l)2C2 


-  i  . 


For"more  anisotropic"  materials  with  m^  <  m  <  m^  (m^  replaced  by  0  if 
m^  <  0),  the  parameters  a^and  B^are  given  by 


b(i)  /  " 


cos  e(i)+  sin  e(i)) . 


with 


I  1  "  K  B(i)  |  > 


4  tan 


'  rnr—  -1 


The  solution  of  the  equations  of  motion  given  by  equations  [6,  and 
[25]  may  be  applied  directly  to  each  layer  of  the  cylinder.  Thus,  the 
complete  solution  for  the  three-layered  cylinder  will  contain  twelve 
integration  constants  A  B  ^  {  (i*l,2,3;  j«*l,2).  These  constants  may  he 

<]  J 

evaluated  by  requiring  the  solution  to  satisfy  the  following  boundary  and 
interface  conditions  of  the  cylinder. 


0  at  r 


0  at  r 


(j-r  ,z) 
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at  r  -  b  if  1-1, 


trJ(1)  -  trJ<1+1)  or  t  -  c  If  1.  2, 


[33  cont] 


(j-r,*). 

By  substitution  of  the  components  of  displacement  and  stress  into 
equations  [33],  a  set  of  twelve  homogeneous,  linear,  algebraic  equations  are 
obtained.  For  a  non-trivial  solution,  the  determinant  of  the  coefficients 
of  the  A^^and  B  ^  must  vanish  resulting  in  the  following  frequency 
equation 


|C  -  0 

ij 


(i.J-1,2 . 12). 


The  non-zero  elements  of  the  determinant  are  given  as 


,  a  (1)  (1)  2  (1)  -  _  .  (IK 

2d  e1  ;  a  a  Z  (;  jj.  a)  - 


r3  w*1*!2  7  2 

;  dA  *  fl(l)  (1)  (l)2  ,  .(1)  ‘>2 

*717(7)  <di  ei  s.  ♦  d5  — ry 

dU)+  dU>  i  5  (1)2  ^ 

5  4  a 


(1)  (*)- 
-dp)Zo(ca  a), 


2  d c2  o(1)  7  U  (;oa)  7)  - 

b  -  1 


.3  .(D-2 

C  dA  a  , .(1)  (1)  (l)2  .(1)  ft2 

da>+  „<i>  <di  ‘i  2  “5  >*7 

J  .  4  ®  S 


.<lh  u  tr  dK 
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j  o 
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C14-2d61)<d51)+<1))^W1  «£(1)a>- 

,U>2  "2  ,.(1)  (DJI)2  (1)  (1)_ 

•d4  «  a  (di  e!  1  +  -T-T-,  -  dc  >  W„  •). 


,(ir  c 


2  .2  5  o 


(1)  2 — 

«  '  •  7^m  •  S?  -  'S’X'V-y. 


df1*  fl2 


•5  '  “4 


(l)  ; 


,(i)  2—  (i)  2 

-4,1^  (d(i),a».w2.  ^4_JL 


„<»♦  d™  1  1  * 

5  4 


(l)2  c2 


-  d^1));a(1)a  W1(ca(1)a), 


2  dfX)  ft2 


Co,  “  "  dl1);a(d.(1)e'[1)6(1)  -  Jt 


12“ 


.(!)*  c‘ 


-  d^)c^!(1)I  zi(ci(1)I), 


9  A(1>  o2 

c24  -:i)u<d?)«l1)l<1)  -  ^T7-^1,)^w1Ui(1)I), 


Cj^»  C4  A  (i  ■  1,2, 3, 4)  -  ^  (i  »  1,2, 3, 4)  with  a  replaced  by  b. 


C3ti*  C4  (i  ■  :.,6,7,8)  ■  C2  i#  ^  (i  ■  1, 2,4,5)  with  a  replaced  by  b 


and  (1)  by  (2), 


C51  -  -  e^;2  a(1)b  Zi  (;a(1)b)  , 


[ 35  cont ] 


r  -  2  JDr  u  ,,  (l)r, 

Cj2  C  5.  0  (;a  b) , 


(1)  (1)  (1)  (1)_ 

C53  -  -  (d^i;  +  d<1})  e^V  Z^B  b). 
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-  -  (d^X)+  d<1})  ;S(1)b  Wx  (;8(1)b). 


5  4 

;2  j—  ,,(1)  (i)  (D2  dr« 


/ 1  \  /n  (d  c  a 

d<i)+  da)  i  i  - 

5  4 


(D.2 

2  2-*d<1))  Z  Ua(1)b), 


,(D  4 
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(d^M”  2(1>2  W  (;a(1)b). 


d(l).  .(1)  '“l  1  - 

d5  +d4 


a 


(i)2  ;2  4  0 


'63 


'64 


,b(d^c^  a(1)2-^4+df>)  zo  (c8(1)b), 
1  2  Cl)2  42  4  0 


M)  2 

—  (1)  (1)  (l)2  d4  W  (1)  H)_ 

4b  (d  C  eu;  -  —  22+  d4  )  Wo  <<1U  »>). 

2  2  *  4 


[ 35  cont ] 


C_  d*  d  (J  "  5.6.7,8)  -  C  ,  C  (i  -  1,2, 3, 4)  with  (1)  replaced  by  (2), 

5,J  6,J  3,1  otl 

C7  d  •  Cft  4  ^  “  5,6, 7,8)  ■  c6  ,,  C  (i  -  1,2, 3, 4)  with  b  replaced  by  c 
/,J  *»,J  “i1  3,1 

and  (1)  replaced  by  (2), 

c7  ,,  CB  4  0  ’  9,10,11,12)  •*  C,  ,  C  (i  -  1,2, 3,4)  with  b  replaced  by 
'•J  8»J  o,i  5, i 

c  and  (1)  replaced  by  (3). 

Co  d’  Ci n  d  ^  “  5,6, 7,8)  "  c.  .  »  c-  .  (j  “  1, 2,3,4)  with  a  replaced  by 

*»J  lu»j  1,J  2,j 

c  and  (1)  by  (2). 

Cg  , »  C1Q  (j  ■  9,10,11,12)  -  C  ,  C  (J  ■  1,2, 3,4)  with  a  replaced 
y  »J  *U*J  l.j  2,j 


by  c  and  (1)  by  (3). 
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CU,J’  C12  J  <J  "  9»10’U'12)  “  C2,j  Ci,j  (J  “  1.2.2.4>  with  a 

_  [ 35  cont ] 

replaced  by  d  and  (1)  by  (3). 

For  given  material  properties  and  shell  geometry,  the  frequency 
equation  [ 34 J  iJ  a  transcendental  relation  between  the  non-dimensional i zed 
frequency  ft  and  the  wave  number  For  any  value  of  the  frequency 
equation  will  yield  an  infinite  number  of  values  of  ft,  each 
corresponding  to  a  different  mode  of  wave  propagation. 

The  frequency  equation  may  be  specialized  to  give  the  correct  formu¬ 
lation  for  solid  rods  (a  ■  0).  In  this  case,  the  boundary  conditions  at 
the  inner  surface  (See  Eq.(33])are  omitted.  Furthermore,  in  order  for 
the  displacement  field  in  layer  (1)  to  remain  finite  at  r  ■  0,  it  must 
contain  no  Bessel  functions  of  the  second  kind.  This  implies  that  the 
constants  of  integration  ^  ( j  ■  1,2)  must  vanish.  Consequently,  the 
frequency  equation  for  axisymmetric  waves  in  three-layered,  transversely 
isotropic  rods  may  be  obtained  by  deleting  the  first  and  second  rows  and 
the  second  and  fourth  columns  of  equation  [34]. 

For  "more  anisotropic'*'  materials,  within  the  range  of  values  of  w 
and  4  for  which  the  values  of  £  and  £  are  complex,  the  frequency  equation 
will  contain  complex  elements.  Inasmuch  as  £  is  the  negative  complex 
conjugate  of  £,  the  following  relations  are  valid 

[  3b  j 

n(ar)*  <-l)n  Y*  (£  r), 

*  * 

where  J  and  Y  are  the  complex  conjugates  of  .1  and  Y  . 


(£  r)  -  (-1 )n 

*  /  — V 

J„  <E  r>> 

<£  7)  «=  (-1)” 

v‘  (£  O. 

n 


n 


n 
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If  the  layer  of  the  cylinder  is"more  anisotropic"  and  if  ra  is 

sufficiently  small  so  that  £  and  6,  are  given  by  equation  [19],  then  the 

th 

four  columns  of  Eq.  [34]  relating  to  the  i  layer  are  complex.  It  can 
be  shown,  that  with  the  exception  of  a  constant  multiplication  factor, 
columns  (4i-n)  and  (4i-n+2)  (n-2,3)  are,  element  by  element,  complex 
conjugate  pairs.  For  instance,  if  layer  1  is  complex  (isl),  the  frequency 
equation  may  be  written  as 


C.  -  G  G  G  G  C  -0 
1  Jk1  1  jl  J2  J3  j4  > 


[37] 


where  C  represents  the  last  eight  columns  of  the  determinant  and  G^ 
(k»l,2,3,4)  are  the  complex  first  four  columns  of  C^.  These  may  be  written 


as 


°J1  ■  V  +  1  V  CJ3  ’  V  ‘ 


J3 

(J-1,2 . 12)  [38] 


G  -  C  +  i  (T  , 
J2  J2  j4 


GJ*  '  -  1  V 


Substitution  into  equation  [37]  and  expansion  yields 


ZA.  \  -  -4  I  Z  t  ?  C  C 
jk1  1  jl  j2  J3  j  4 


-  0 


[39] 


Thus,  the  problem  of  evaluating  the  frequency  equation  for  waves  traveling 
in  shells  made  of  "more  anisotropic"  materials  is  essentially  identical  to 
that  of  evaluating  the  frequency  equations  for  waves  traveling  in  shells 
made  of  "less  anisotropic"  materials. 
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WAVES  WITH  INFINITELY  LON-.  AXIAL  WAVE  LENGTHS 


When  the  axial  wave  number  vanishes,  the  displacement  field  of  the 
cylinder  is  independent  of  the  axial  coordinate.  The  radial  wave  numbers 
reduce  to 


,(i) 


sr 


.  (i)‘ 


and 


2  d(i> 
(i)2  d4 


TiT  p 


(i) 

c 


[40] 


a  '  1 

and  the  frequency  determinant  can  be  written  as  the  product  of  the 
following  two  determinants 


where 


D  - 


and 


V 


Limit 

c-*o 

'V 

d2-o 

• 

C21 

C22 

0 

0 

0 

0 

C31 

C32 

C35 

C36 

0 

0 

C61 

C62 

C65 

C66 

0 

0 

0 

0 

C75 

C76 

C79 

n 

h- 

o 

0 

0 

C10,5 

C10,6 

C10, 9 

C10,l( 

0 

0 

0 

0 

C 

11,9 

c 

11, 1( 

C13 

C14 

0 

0 

0 

0 

C43 

C44 

°47 

C48 

0 

0 

C53 

C54 

C57 

C58 

0 

0 

0 

0 

C87 

C88 

C8,ll 

C8,12 

0 

0 

C97 

C98 

C9,ll 

C9,12 

[41] 


[42] 


[43. 


0 


0 


0 


0 


C12, 11  C12, 12 
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The  equation  *  0  yields  the  cut-off  frequencies  of  axisymraetric 
longitudinal  shear  vibrations  involving  only  axial  displacement.  In  this 
case,  the  motion  is  uncoupled  equivoluminal,  and  the  displacement  components 


u(i)  -  0, 
r  * 


“z1’  ■  -»ci)  l  4l)  >0  (»=i>r)  +  “l1’  Y0  (pc1>r) 


sin  ft  t. 


Notice  that,  as  in  the  case  of  isotropic  shells,  the  motion  depends 
only  on  the  elastic  constants  d^(i=l,  i,  3). 

The  equation  D  .  =  0  represents  plane  strain  extensional  motion  involving 
only  radial  displacements 

u*i}  -  q^  [A^  Jx  (q^r)  +  B2  ^  Y1  cos  ft  t. 

r  c^xc  zic  [45] 

u(i)-  0. 
z 

This  motion  is  independent  of  the  elastic  constants  d^^  and  d^  and 
consequently  is  independent  of  the  axial  Young's  moduli  of  the  three 


materials. 


i 

t 


NUMERICAL  ANALYSIS 


A  computer  program  has  been  written  for  numerical  evaluation  of  the 
frequency  equations.  The  program  first  computes  the  cut-off  frequencies 
on  the  basis  of  Eqs.  [42]  and  [43],  and  utilizes  them  as  starting  values 
to  trace  the  branch  curves  of  Eq.  [34]  on  the  ft  -  c  plane.  For  each 
assumed  value  of  c,  the  frequency,  ft,  is  incremented  by  a  specified  amount 
Aft  until  a  change  in  the  sign  of  the  determinant  occurs.  This  indicates 
a  root  between  the  last  two  values  of  ft.  An  interval  halfing  procedure 
is  then  executed  which  brackets  the  root  to  the  required  (pre-set)  accuracy. 
Subsequently,  the  value  of  c  is  incremented  by  a  pre-assigned  increment 
Ac  and,  starting  with  a  new  value  of  ft,  (computed  from  the  slope  of  the 
two  previously  established  points  on  the  branch)  the  process  is  repeated 
and  new  roots  are  established  until  each  branch  of  the  frequency  equation 
has  been  traced  up  to  a  pre-assigned  value  of  c.  For  each  tested  value 
of  ft  and  c,  the  program  uses  the  appropriate  form  of  the  frequency  equation 
depending  on  whether  the  radial  wave  numbers  p  and  q  are  real,  imaginary 
or  complex. 

(2)  (2) 

The  effect  of  the  ratio  of  elastic  constants  c. ,  and  c, ,  (for 

4  ♦  11 

constant  a  ratio  of  c^^/  c^^  ■  7/3)  on  the  cut-off  frequencies  of  Lhe 
first  few  modes  is  illustrated  in  Figs.  3  through  6.  In  these*  figures,  the 
outside  and  the  inside  layer  were  chosen  as  Boron/Epoxy  (see  Table  1). 
This  material  is  "less  anisotropic".  The  results  shown  are  valid  for  any 
values  of  ^  and  Cj^^  inasmuch  as  the  cut-off  frequenciesf equations 

[42]  and  [43]]are  independent  of  these  elastic  constants.  As  expected,  the 
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frequencies  of  the  longitudinal  shear  modes  are  independent  of  the  elastic 
(2)  2 

constants  c^j  and  Cj^  .  Moreover,  as  can  be  seen  from  Fig.  3  and  4, 
the  frequency  of  the  first  longitudinal  shear  mode  is  not  noticeably  effected 
by  changes  in  the  axial  shear  modulus  c^^.  The  frequency  for  this 
mode  is  slightly  larger  for  rods  (Fig.  4)  than  for  thick-walled  shells,  i.e. 

H / R  =  1 . 0  (Fig.  3).  For  values  of  H/ R  <  1 . 0  the  effect  of  H/ R  on  the  frequency 
of  this  mode  is  negligible.  Thus,  for  thin  sandwich  shells  the  frequency  of 
the  first  longitudina?  shear  mode  is  only  effected  by  changes  in  the  density 
ratio  of  the  layers.  This  result  is  interesting  inasmuch  as  this  frequency 
is  employed  in  establishing  the  correction  factor  in  Timoshenko-type  shell 
theories . 


TABLE  1.  -  MATERIAL  CONSTANTS 


Material 

6 

Pounds  Per  Square  Inch  (10) 

3 

p (lb/ in  ’ ) 

C11 

C12 

c 

33 

C44 

C13 

Aluminum 

13.46 

5.76 

13.46 

3.85 

5.76 

0.100 

Boron/ Epoxy 

3.28 

1.19 

30.4 

1.00 

0.995 

0.075 

Composite 

Beryllium 

4.24 

0.388 

4.88 

2.36 

0.203 

0.067 
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Figure  7  shows  the  frequency  spectrum  of  a  sandwich  shell  made  of 
an  aluminum  core  and  fiber-reinfor ced  composite  facings  made  of  an  epoxy 
matrix  reinforced  by  unidirectional  boron  fibers.  Figure  8  shows  the 
frequency  spectrum  of  a  sandwich  shell  made  of  an  aluminum  core  and 
beryllium  plate  facings.  The  mechanical  properties  of  the  layers  of  these 
shells  are  given  in  Table  1.  The  aluminum  is  isotropic,  whereas  the 
Boron/Epoxy  composite  and  the  beryllium  are  less  "anisotropic"  and 
"more  anisotropic"  respectively.  The  ft-£  plane  may  be  sub-divided  into 

ml1) 


sectors  by  lines  ft  =  b^£,  ft  =  a^£,  and  ft  = 


Hi=  1.2,3;  j  =  1,2). 
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Throughout  each  sector  the  radial  wave  numbers  p  and  q  retain  their  real, 

imaginary  or  complex  character  and,  consequently,  the  form  of  the  frequency 

equation,  in  each  sector,  does  not  change.  It  can  be  seen  that  the  behavior 

of  the  frequency  lines  of  the  lowest  two  modes  of  wave  propagation  in 

the  shell  having  fiber-reinforced  composite  facings  differs  considerably 

from  those  for  the  shell  with  the  beryllium  facings;  those  of  the  latter 

cross  into  the  region  of  complex  p  and  q.  The  frequency  lines  of  the 

higher  modes  for  the  two  shells  are  comparable  for  large  wave  numbers;  they 

appear  to  become  parallel  to  the  ft  ■  a v  \  line.  For  smaller  wave  numbers, 

within  a  certain  range  of  i, ,  the  frequency  lines  for  the  shell  with  fiber- 

(b) 

reinforced  composite  facings  become  nearly  parallel  to  the  ft  *  b  ;  line. 
This  tendency  is  not  apparent  in  the  spectrum  for  the  shell  with  the 


beryllium  facings. 
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it.  abstract 


The  propagation  of  Infinite  trains  of  symmetric  harmonic  waves 
traveling  in  infinitely  long,  right  circular  cylindrical  shells  Is 
Investigated  on  the  basis  of  the  three-dimensional  theory  of  elasticity. 
The  shells  are  assumed  i..ade  of  three  concentric,  transversely  isotropic 
cylinders  each  of  different  materials  bonded  perfectly  at  their 
Interfaces.  The  frequency  equation  is  established  by  representing 
the  displacement  field  in  each  cylinder  in  terms  of  potential  functions 
and  satisfying  the  Navier  equations  of  motion  and  the  boundary  and 
interface  conditions  of  the  cylinder. 

The  frequency  equation  has  been  programmed  for  numerical  evaluation 
on  an  IBM  70^4/7094  DCS  computer,  and  the  influence  of  the  mechanicel 
properties  of  the  layers  on  the  frequencies  of  the  first  few  modes  is 
i  nvest  l  gated. 
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